Quasi-Scarred Resonances in a Spiral-Shaped Microcavity 
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We study resonance patterns of a spiral-shaped dielectric microcavity with chaotic ray dynamics. 
Many resonance patterns of this microcavity, with refractive indices n = 2 and 3, exhibit strong 
localization of simple geometric shape, and we call them quasi-scarred resonances in the sense 
that there is, unlike the conventional scarring, no underlying periodic orbits. It is shown that 
the formation of quasi-scarred pattern can be understood in terms of ray dynamical probability 
distributions and wave properties like uncertainty and interference. 
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The scar phenomenon, since its advent in a chaotic 
biUiard, has attracted much attention because it had 
not been anticipated from the prevailed random matrix 
theory 0. It is now known that the scarred eigenfunc- 
tions show not only strong enhancement along a unstable 
periodic orbit, but also detail of the stable and unstable 
manifolds around the periodic orbit 3]. This scar effect 
therefore has been regarded as an important feature of 
chaotic systems different from random systems. Another 
important aspect of the scarring effect is its ubiquitous 
existence; it has been observed in various chaotic sys- 
tems such as microwave cavity 0], semiconductor quan- 
tum well 0, surface wave 0, optical cavities etc. 

Recently there is a considerable interest in the light 
emission from dielectric cavities with chaotic ray dynam- 
ics, since many intriguing light emission behaviors take 
place and are known to be relevant to the underlying 
chaotic ray dynamics . There are several reports of ob- 
servation of scarred lasing modes in dielectric microcav- 
ities of various boundary shapes 0, IE 0| ■ The scarred 
lasing modes generally show good directionality of light 
emission, the directionality is an important characteris- 
tic required for applications to photonic and optoelectric 
information processing [h^. The number of directional 
beams of the scarred emission from usual microcavities 
would be more than two because of the discrete sym- 
metry of the boundary geometry and the possibility of 
interchanging incident and reflected rays on the underly- 
ing periodic orbit. 

In a remarkable experiment, Chern et al. have success- 
fully observed unidirectional emission in spiral-shaped 
quantum- well microlasers [T^ . The unidirectional laser 
beam is important to arrange easy optical communica- 
tion between microlasers. The spiral-shaped boundary, 
in which ray dynamics is chaotic, is given by 



r(0) = R{1 + 



(1) 



in polar coordinates (r, 0), where R is the radius of the 
spiral at = and e is the deformation parameter. Ba- 
sically, the unidirectionality of the emission beam comes 
from the special properties of the spiral-shaped bound- 
ary geometry which other common cavity designs do not 
have. They are the absence of any symmetry and the ex- 



istence of the notch. As mentioned above, the absence of 
symmetry would be the necessary condition for the uni- 
directional emission. The notch makes the microcavity 
show very strong chirality by transmitting or reflecting 
counterclockwise rotating rays. The bouncing from the 
notch is inevitable for the rays and would be an essen- 
tial process for the unidirectional emission. Besides the 
unidirectionality, it is important and interesting to study 
how the unique characteristics of the spiral-shaped mi- 
crocavity appear on resonance patterns. 

In this Letter, we investigate the resonance patterns in 
the spiral-shaped dielectric microcavity. We find that a 
large number of resonances obtained are strongly local- 
ized and that the localized patterns are not supported 
by any unstable periodic orbit, so we call them quasi- 
scarred resonances. The existence of quasi-scarred reso- 
nances implies that the scarring phenomenon in dielectric 
microcavities has substantial differences from the conven- 
tional scarring in billiard systems. The differences come 
from inherent characteristics of dielectric cavities such 
as existence of the critical incident angle for total in- 
ternal reflection and energy loss by refractive emission. 
We explain the formation of the quasi-scarred resonances 
in terms of ray dynamical probability distributions and 
wave properties like uncertainty and interference. For 
convenience, we take e = 0.1 and i? = 1 in this Letter. 

In order to investigate the ray dynamical properties 
of the spiral-shaped dielectric microcavity, we first con- 
sider a uniform ensemble of initial points over the whole 
phase space (s, p), where s is the boundary arc length 
from the = point (see Fig. 4) and its conjugate 
variable p is given as p = sinO, 9 being the incident an- 
gle of ray. If the boundary is made by a perfect mirror, 
the distribution of the points in the phase space at later 
times would remain uniform (in a random sense) and 
structureless. However, in the dielectric microcavity, the 
distribution of the points is, some time later, not uni- 
form but rather structural because the individual ray can 
suffer energy loss by refractive emission when bouncing 
from boundary. The amount of the energy loss is de- 
termined by the transmission coefficient T{p) 13] which 
has a nonzero value in the range of ~pc < p < Pc, where 
Pc is the critical line for total internal reflection and is 
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FIG. 1: The escape time distributions Pes(t) with n = 2 for 
two different sets of initial points; the thick and thin hnes 
represent numerical results for Set A and Set B, respectively. 
The time t is scaled to be the length of ray trajectory when 
R=l. 

related to the refractive index n as pc = sin 6'c = 1/n, 
9c being the corresponding critical incident angle. This 
leaky property of rays in the ensemble is described by the 
survival probability distribution P(s,p, t), the probability 
with which the ray with (s, p) can survive in the micro- 
cavity at a time t. With the survival probability distribu- 
tion P{s,p, t), the energy £{t) confined in the microcavity 
and the escape time distribution Pes{t) are expressed as 
£{t) — £o J ds dp P{s,p,t), £q being the initial energy, 
and Pes(t) — j ds dp P{s,p,t)T{p), respectively. Since 
the confined energy decreases by the ray transmission 
through cavity boundary, we can get a relation. 



d£{t) 
dt 



-SoPesit). 



(2) 



It is well known that in fully chaotic open systems 
the escape time distribution Pes{t) shows exponential 
long time behavior, while it becomes power law decay 
in the KAM systems due to the stickiness of the KAM 
tori The exponential behavior of Pes{t) suggests 

that P{s,p, t) would have the same phase space distribu- 
tion after a certain period of time, i.e.. 



Pis,p,t) = B{t)P,{s,p), 



(3) 



which defines the steady probability distribution Ps{s,p) 
as the stationary part of P{s,p,t). It is obvious from 
Eq. |5J) that the relation in Eq. (j^J is equivalent to as- 
suming the exponential time behaviors of ray dynamical 
distributions such as £{t), P{s,p,t), and Pes{t)- In the 
case of dielectric microcavities, a numerical justification 
of the relation in Eq. ^ will be presented below (see 
Fig. 1). The steady probability distribution Ps{s,p) then 
characterizes the ray dynamical long time behavior. The 
decay rate 7, from Eq. can be expressed as 



7 



dsdpPs{s,p)T{p), 



(4) 




and the ray dynamical near field and far field distribu- 
tions can be also described by Ps{s,p). 
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FIG. 2: (a) The normalized survival probability distribution 
P{s,p, t) sampled in the time range of 57 < t < 60 for Set A. 
The white dots represent points whose weight is greater than 
0.1 after normalization, (b) The distance, resulting after 3 
bounces, distribution d{s,p). From comparing both figures, 
it is clear that the partially reflected rays near —pc would 
make rough triangular trajectories. 

For simplicity's sake, we will concentrate on TM 
(transverse magnetic) polarization in this Letter. In 
Fig. 1 , the escape time distributions Pes (t) are shown for 
the n = 2 case. Here, we consider two different sets of 
400 X 400 initial points: one is the uniformly distributed 
set over the whole phase space (Set A) and the other is 
the uniformly distributed one in a part of the phase space, 
(0 < s < Smaxfi, 0.5 < p < 0.75) (Set B), where Smax 
is the total length of the boundary. Note that above 
tc ~ 30 exponential decay behaviors are shown. The 
slope of the linear part determines the decay rate 7. The 
similar slopes for both Set A and Set B reflect that rays 
lose their energy through the same process. The details 
of the process appear in the structure of Ps{s,p). 

Figure 2 (a) shows an approximate Pg(s,p) for n ~ 2 
given by normalizing the P{s,p,t) in the time range of 
57 < i < 60 for Set A. The structure of the approximate 
Ps (s, p) is almost invariant in other time ranges of the lin- 
ear part {t > tc) and even for Set B. It is clear that the 
energy loss was mainly caused by tangential emissions 
just above the critical line {—pc = ~l/2)- So, we can 
see that the process mentioned above is the way that the 
ray trajectories first rotate counterclockwise {p > Pc), 
then change their rotational direction by reflection on 
the notch part, and afterwards gradually approach ~Pc- 
Most of them are emitted out from the microcavity and 
the remains repeat the same process. The distribution 
confined to the negative value of p means strong chiral- 
ity of this spiral-shaped microcavity. The dark tentacu- 
lar structure in Fig. 2 (a) implies the missing trajectories 
which are reflected at the notch with \p\ > pc- In fact, 
the overall structure presents a part of unstable mani- 
folds, which is typical in open chaotic systems This 
structure would give important informations about sta- 
tistical properties of resonances, i.e., far field and near 
field distribution of resonances would show minima at 
values corresponding to the missing trajectories. 

More direct implication on resonance patterns 
can arise from the distribution of resulting dis- 
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FIG. 3: (color). Field intensity plots of quasi-scarred reso- 
nances in the spiral-shaped microcavity. (a) n = 2 and nkR = 
(109.70,-0.1128). (b) n = 3 and nkR = (109.59,-0.1127). 
In figures, the field intensity is normalized by scaling the max- 
imum intensity as one. 

tance after 3 bounces (Fig. 2(b)), i.e., d{s,p) — 
\/{sf ~ s)2 + {pf ~ p)"^ where (s, p) is the initial posi- 
tion and (s/, Pf) being the position after 3 bounces. We 
note that in Fig. 2 (b) the critical line p — —pc lies on 
the region of lower d values. Since the rays in the region 
just above —pc are partially emitted out, the remaining 
reflected rays would make a rough triangle. As discussed 
below, the imprint of this fact appears apparently in res- 
onance patterns (see Fig. 3 (a)). Although the n = 3 case 
is not presented in the figures, the distance distribution 
after 5 bounces also shows similar features, implying that 
the star shape ray trajectories would be responsible for 
resonance patterns (see Fig. 3 (b)). 

Using the boundary element method p^. we obtain 
resonances around Re{nkR) ~ 110 for the spiral-shaped 
dielectric microcavity, 24 resonances for n = 2 and 23 
resonances for n = 3, which are about 25% of the to- 
tal number of resonances in the concerned range. From 
the resonances we realize an important fact that the ba- 
sic localized structures of the resonance patterns are tri- 
angular and star shapes for n — 2 and 3, respectively, 
which is consistent with the implication of Ps {s,p). Some 
of these resonances look like strongly scarred eigenfunc- 
tions of billiard system, showing strong directional emis- 
sions matched to the triangular and star patterns. The 
nkR values and patterns for whole resonances will be 
presented elsewhere due to lack of space. 

The most clearly localized resonances for n = 2 and 
3 are shown in Fig. 3. The patterns look like strongly 
scarred resonances, but there is no exact underlying un- 
stable periodic orbit. Absence of periodic orbits of simple 
geometry, without bouncing at notch, e.g., triangle and 
star, is evident by numerical evaluation of dp = pi—pf for 
a closed triangle or star trajectory starting from {si,pi) 
and terminating at (s^, pf). Moreover nonexistence of 
periodic orbits of simple geometry can be understood if 
one know that for the clockwise rotating case the distance 
between origin and the ray segment always decreases as 
far as the ray bounces at the curved part of the spiral- 



FIG. 4: Schematic diagram for quantifying the degree of 
uncertainty. The trajectory satisfying the Snell's law, with 
an incident angle (012 + (/>i3)/2, is denoted by dashed lines. 

shape boundary. We obtain \5p\ > a for arbitrary Si 
value, where a = 0.075 for the triangle trajectory and 
cr = 0.136 for the star trajectory. Since the localized pat- 
terns of resonances are not supported by any unstable pe- 
riodic orbit, we call them quasi-scarred resonances. The 
existence of quasi-scarred resonances in dielectric cav- 
ities can be understood from the inherent property of 
dissipative systems, i.e., uncertainty characteristics. An- 
other important result from the resonance pattern analy- 
sis is that many resonances are quasi-scarred, e.g., in the 
present case more than a half are quasi-scarred, while 
only a small fraction of eigenfunctions are scarred in bil- 
liard systems. In practical experiment, this implies that 
the quasi-scarred lasing emission can be excited easily 
due to its dominant existence in resonances. In fact, 
the dominant existence of quasi-scarred resonances can 
be regarded as a result of the openness of microcavities. 
In open systems, rather local part of phase space would 
support resonances(e.g., see Fig. 2(a)) and the result- 
ing individual resonance would show a stong localization 
whose pattern might be determined by the property of 
the openness. This is consistent with results, associated 
with scarred resonances, in various open systems 0, 

Now, we consider bouncing positions of the triangle 
formed in quasi-scarred resonances which seem to have 
a definite dependence on their Re(nfci?) values. We as- 
sume that the triangle in quasi-scarred resonances has 
minimum deviation from the ray trajectory governed by 
the Snell's law, and maximum constructive interference 
under constraint of high intensity of the electric field at 
the bouncing positions. We quantify these by two factors, 
a and (3 as follows. Let Si [i — I, 2, ?>) be the bounc- 
ing positions of a triangle, from the angles {4>ij, 4>ik) 
to the normal line on the boundary, and we can define 
Pij = sm{(j)ij),pik = sin(0ifc), andpi = sin(^i^^^) (here 
i, j, k are cyclic). Also we get the new positions Sij, Sik 
as the next positions of (si, pi) and (si, ~ Pi), respec- 
tively (see Fig. 4). Then we define partial uncertainty of 
the triangle given by (si, S2, S3) as 

Oil = [[pi - Ptj){sj - s.ij)f + [{pi - Ptk){sk - Sik)Y. (5) 

Total uncertainty, therefore, is the sum of these terms, 
a — 'Yl^i=i^i- definition, when the triangle is a 

periodic orbit, a becomes zero. To quantify the de- 
gree of constructive interference we consider mi -\- (3i = 
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FIG. 5: Variation of the optimized bouncing positions 
(si, S2, S3). The solid lines denote the present theory with a 
correction /i = 0.06. The circles represent the bouncing posi- 
tions of triangular quasi-scarred resonance patterns of n = 2 
case; three solid circles with the same Ke{nkR) correspond to 
the main triangular pattern, and three open circles do to the 
secondary triangular pattern in a quasi-scarred resonance. 

li/{X/2) + S(f)/TT for each triangle segment of length k, 
where rrii is an integer and —0.5 < [3i < 0.5, A — 27r/(nfc), 
and 54> is the phase shift arisen from total internal reflec- 
tion jl 7j . Total quantity for the degree of constructive 
interference is then (3 — 0i with an additional con- 

straint that the sum M = X]i=i should be even. 

We first determine triangles with minimum uncertainty 
a as a function of si, and then apply the condition 
of minimum [3 to the triangles. From this process we 
get the most optimized triangle of {s\^ S2, S3) for a 
fixed Re(nfci?). The direct application of this method 
shows systematic deviation from bouncing positions of 
resonance patterns. This systematic discrepancy results 
from the fact that rays inside microcavity have angu- 
lar distributions and, also the boundary has curvatures, 
which gives rise to a correction of the Snell's law. This 
effect is prominent near the critical angle 9c^ studied and 



known as Goos-Hanchen and Fresnel Filtering ef- 
fects 0, 01 • We here incorporate these effects by taking 
effective segment length I* = U + fiX. The results are 
shown in Fig. 5. The solid lines are results of the present 
theory with fi = 0.06 which are in good agreement with 
the bouncing positions (denoted by circles) of the quasi- 
scarred resonances. Absence of bouncing positions near 
s = 2.0 and s = 4.5 is consistent with the tentacular 
structure of the approximate steady probability distribu- 
tion Ps{s,p) in Fig. 2 (a). 

In conclusion, we have found that the localized pat- 
terns of resonances in a spiral-shaped dielectric microcav- 
ity are constructed by the quasi-scar phenomenon which 
comes from inherent properties of the dielectric micro- 
cavity, and that a large fraction of the resonances are 
quasi-scarred. The results are contrasted with the case 
of billiard systems in which only scar phenomenon exists, 
and a small fraction of eigenfunctions are scarred. Even 
though the system is chaotic, it is possible to extract 
some information on resonance patterns from the ray dy- 
namical consideration, more precisely, from the steady 
probability distribution Ps{s,p). Since Ps{s,p) contains 
long lasting ray dynamical information, its structure 
should be related to the high-Q resonances which are 
likely to appear as lasing modes. From a theoretical 
viewpoint, just like the semiclassical approach in Hamil- 
tonian systems 20|, a semiclassical method in dielectric 
cavities might be useful to understand resonance posi- 
tions and degree of scarring or quasi-scarring. Develop- 
ing semiclassical theory for dielectric cavities seems to 
be nontrivial. We expect that the results of this Letter 
will improve physical insight onto resonance patterns in 
microcavities. 
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